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ABSTRACT 

We perform evolutionary calculations of binary stars to find progenitors of systems 
with parameters similar to the recurrent novae U Sco. We show that a U Sco type- 
system may be formed starting with an initial binary system which has a low-mass 
carbon-oxygen white dwarf as an accretor. Since the evolutionary stage of the sec- 
ondary is not well known, we calculate sequences with hydrogen rich and helium 
rich secondaries. The evolution of the binary may be devided into several observable 
stages as: classical nova, supersoft X-ray source with hydrogen stable burning and 
strong wind phases, ending up with the formation of a massive white dwarf near the 
Chandrasekhar mass limit. We follow the chemical evolution of the secondary as well 
as of the matter lost from the system, and we show that observed 12 C/ 13 C and N/C 
ratios may give some information about the nature of the binary. 



Key words: 

lution 



binaries: close 
stars: recurrent novae 



binaries: general stars: mass loss evo- 

star: individual: U Sco 



1 INTRODUCTION 

Recurrent novae are a small class of objects which bear many 
similarities to other cataclysmic variable systems. They ex- 
perience recurrent outbursts at intervals of 20-80 yrs. 

Webbink et al. (1987) lengthily discussed the nature 
of the recurrent novae, and they concluded that according 
to outburst mechanisms there are two subclasses of these 
systems: (a) powered by thermonuclear runaway on the sur- 
face of the white dwarf (e.g. U Sco), and (b) powered by 
the transfer of a burst of matter from the red giant to the 
main-sequence companion. 

U Sco is one of the best observed recurrent novae. His- 
torically, its outbursts were observed in 1863, 1906, 1936, 
1979, 1987 and in 1999. Determinations of the system vi- 
sual luminosity at maximum and minimum, indicate a range 
Amy ~9. We also note that the mean recurrent interval im- 
plied by known outbursts is P rec = 23 yrs. Schaefer (1990) 
and Schaefer & Ringwald (1995) observed eclipses of U Sco 
in the quiescent phase, and determined the orbital period 
P ori) =1.23056 d. 

Ejecta abundances have been estimated (from 1979 out- 
burst) from optical and UV studies by Williams et al. (1981) 
and Barlow et al. (1981). They derived extremely helium rich 
ejecta He/H~ 2 (by number), while the CNO abundance 
was solar with an enhanced N/C ratio. From the analysis of 
the 1999 outburst Anupama & Dewangan (2000) obtained 



an average helium abundance of He/H~0.4±0.06. The esti- 
mated mass of the ejected shell for 1979 and 1999 outbursts 
is ~ 10~ 7 Mq (Williams et al. 1981; Anupama & Dewan- 
gan 2000). Spectroscopically, U Sco shows very high ejection 
velocities of (7.5-11) x 10 3 km s _1 (Williams et al. 1981; Mu- 
nari et al. 2000). Latest determinations of the spectral type 
of the secondary indicate a K2 subgiant (Anupama & De- 
wangan 2000; Kahabka et al. 1999). Following to Kahabka 
et al. (1999) the distance to U Sco is about 14 kpc. 



According to Kato model (1996), supersoft X-ray emis- 
sion is predicted to be observed about 10-60 days after the 
optical outburst. BeppoSAX detected supersoft X-ray emis- 
sion from U Sco in range 0.2-20 keV just 19-20 days after 
the peak of its optical outburst in February 1999 (Kahabka 
et al. 1999). The fact that U Sco was detected as a super- 
soft X-ray source (SSS) is consistent with steady hydrogen 
burning on the surface of its white dwarf component. 



In this paper we construct several evolutionary se- 
quences which may lead to formation of systems like U Sco. 
In Section 2 the evolutionary code is briefly described. In 
Section 3 we discuss the effect of mass transfer on binary 
evolution. Section 4 contains the results of the calculations. 
A general discussion and conclusion follow. 
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2 THE EVOLUTIONARY CODE 

The models of secondary stars filling their Roche lobes were 
computed using a standard stellar evolution code based on 
the Henyey-type code of Paczyhski (1970), which has been 
adapted to low-mass stars (Marks & Sarna 1998, hereafter 
MS98). 

Our nuclear reaction network is based on that of 
Kudryashov & Ergma (1980), who included the reactions 
of the CNO tri-cycle in their calculations of hydrogen and 
helium burning in the envelope of an accreting neutron star. 
We have included the reactions of the proton-proton (PP) 
chain. Hence we are able to follow the evolution of the el- 
ements: 1 H, 3 He, 4 He, 7 Be, 12 C, 13 C, 13 N, 14 N, 15 N, 14 0, 
15 0, le O, 17 and 17 F. We assume that the abundances of 
18 and 20 Ne stay constant throughout the evolution. We 
use the reaction rates of: Fowler, Caughlan & Zimmerman 
(1967, 1975), Harris at al. (1983), Caughlan et al. (1985), 
Caughlan & Fowler (1988), Bahcall & Ulrich (1988), Bahcall 
& Pinsonneault (1992), Bahcall, Pinsonneault & Wesserburg 
(1995) and Pols et al. (1995). We use the Eggleton (1983) 
formula to calculate the size of the secondary's Roche lobe. 

For radiative transport, we use the opacity tables of 
Iglesias & Rogers (1996). Where they are incomplete, we fill 
the gaps using opacity tables of Huebner et al. (1977). For 
temperatures lower than 6000 K we use the opacities given 
by Alexander & Ferguston (1994) and Alexander (private 
communication). For a more detailed description of the code 
see MS98. 



3 THE EFFECTS OF MASS TRANSFER 

While calculating evolutionary models of binary stars, we 
must take into account mass transfer and associated physi- 
cal mechanisms which lead to mass and angular momentum 
loss. Apart from angular momentum loss due to gravita- 
tional wave radiation and magnetic braking, there is addi- 
tional loss due to the non-conservative nature of semide- 
tached evolution, such as novae outbursts and strong op- 
tically thin/thick wind from white dwarf occuring during 
stable hydrogen shell burning phase. 

We can express the change in the total orbital angular 
momentum (J) of a binary system as 
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where the terms on the right hand side are due to: grav- 
itational wave radiation, magnetic stellar wind braking, no- 
vae outbursts angular momentum loss (which describe the 
loss of angular momentum from the system due to non- 
conservative evolution), frictional angular momentum loss 
(which occurs during a novae outburst as the secondary or- 
bits within the expanding novae shell and during strong op- 
tically thin/thick wind) and optically thin/thick wind from 
white dwarf (during stable hydrogen shell burning). The role 
of the first two terms have been discussed in many papers 
(see, for example, review paper by Verbunt 1993). In our 
paper we discuss the last three terms. 



3.0.1 Angular momentum loss associated with novae 
outbursts 

To take account of the angular momentum loss that accom- 
panies mass loss due to novae outbursts that occur on the 
surface of the white dwarf during the semidetached phase, 
we use a formula based on that used to calculate angular mo- 
mentum loss via a stellar wind (Paczyhski (1967); Ziolkowski 
(1985) and De Greve (1993)), 



M wd Msg 
M 3 g M tot ' 



where M = /1 M sg . (2) 



/1 is the ratio of the mass ejected by the white dwarf 
to that accreted by the white dwarf; /2 is defined as the 
effectiveness of angular momentum loss during mass transfer 
(Sarna & De Greve 1994, 1996); M wd and M sg denote mass 
of white dwarf primary and subgiant secondary, respectively; 
Mtot (=M w d + M S g) is the total mass of the system; M and 
Msg are, respectively, the rate of mass loss from the system 



and the rate of mass loss from the secondary 



Ms 



equivalent to the mass transfer rate). We take / 2 = 1.0 which 
is typical for a stellar wind. Recently, Livio & Pringle (1998) 
proposed a model in which the accreted angular momentum 
is removed from the system during novae outbursts, which 
agrees with our earlier suggestions (Marks, Sarna & Prialnik 
1997, MS98). Equation (2) is in quantitative agreement with 
estimations made by Livio & Pringle (1998). 



3.1 Frictional angular momentum loss 

During nova outburst or optically thin/thick wind phase the 
secondary star effectively orbits within the expanding nova 
shell or dense wind matter. Due to the frictional deposition 
of orbital energy into expanding nova shell or strong wind, 
the separation of the components will be decreased. 

Livio, Govarie & Ritter (1991) estimated the change in 
the orbital angular momentum brought about by frictional 
angular momentum loss over a complete novae cycle, 



1 (1 + U 



2\ 2 



M wir 
M 3C 



(3) 



where U is the ratio of the expansion velocity of the 
envelope at the position of the secondary to the orbital ve- 
locity of the secondary in the primary's frame of reference 
(fcxp/forb), (g = M wd /M 3g ) is the mass ratio and M wind 
is the rate of mass flow past the secondary Warner (1995). 
For strong winds from white dwarf during stable hydrogen 
shell burning Hachisu, Kato & Nomoto (1996) estimated 
^cxp/uorb ~ 10. Since we have no information concerning 
the expansion velocity of the ejecta from theoretical novae 
models of Prialnik & Kovetz (1995) and Kovetz & Prialnik 
(1997) (hereafter PK95 and KP97), we do not include the ef- 
fect of frictional angular momentum loss during nova phase. 
We include this effect during strong wind phase. 



3.2 Optically thin/thick wind angular momentum 
loss 

If we consider a carbon-oxygen (C-O) white dwarf accreting 
matter from a companion with solar composition, there ex- 
ists a critical accretion rate above which the excess material 



© 0000 RAS, MNRAS 000, 000-000 



An evolutionary scenario.. 3 



is blown off by strong wind. Hachisu et al. (1996) show that 
because wind velocity is about 10 times higher than orbital 
velocity, the wind has the same specific angular momentum 
as that of the white dwarf, which is estimated as 



q M wind 
l + q M sg 



(4) 



As argued by MS98, for typical binary parameters, an- 
gular momentum loss due to nova outbursts and magnetic 
stellar wind braking are the dominant angular momentum 
loss mechanisms, with gravitational wave radiation two or- 
ders of magnitude less effective (when the orbital period de- 
creases sufficiently, gravitational wave radiation will be more 
effective). Frictional angular momentum loss is the least ef- 
fective at eight orders of magnitude less than novae out- 
bursts angular momentum loss and magnetic stellar wind 
braking. However this mechanism will be very effective dur- 
ing strong optically thin/thick wind (SSS phase). Note also, 
that during the wind phase the wind carries off the specific 
angular momentum of the white dwarf, which stabilizes the 
mass transfer (Hachisu et al. 1996; Li & van den Heuvel 
1997). 



3.3 Accretion of material ejected during novae 
outbursts 

We employ the code developed by MS98, which utilises the 
results of the theoretical novae calculations made by PK95 
and KP97. By interpolation from the data sets of PK95 and 
KP97, at each time-step, we use the white dwarf mass and 
the mass transfer rate to determine the novae characteristics: 
/i, amplitude of the outburst, recurrence period, chemical 
composition of the ejected material. 

To calculate the re-accretion by the secondary of ma- 
terial ejected during novae outbursts, we assume that the 
mass of the material re-accreted (M re _ acc ) is proportional 
to the mass of the material ejected by the white dwarf such 
that, 



4a 2 



(5) 



outburst, R sg and a are the radius of the secondary star and 
the separation of the system, respectively. The constant of 
proportionality is the ratio of the cross-sectional area of the 
secondary star to the area of a sphere at radius a from the 
white dwarf. We base this formula on the assumption that 
novae ejections are spherically symmetric and instantaneous. 



3.4 Accretion of material from strong wind 

We define two critical mass accretion rate onto white dwarf. 
First (Warner 1995): 



M cr ,i = 2.3 x (M wd -0.19) 3/2 , 



(6) 



describing critical accretion rate above which hydrogen 
rich material is burning in a stable shell; and below which 
novae outbursts occur. Second (Nomoto, Nariai & Sugimoto 
1979; Hachisu et al. 1996): 



describing the critical accretion rate above which strong 
wind solution by Hachisu et al. (1996) is working. It allows 
burning of the hydrogen into helium at a rate close to M C r,2, 
with the excess material being blown off by the wind at the 
rate: 



M wind = M s , 



M c 



(8) 



To calculate re-accretion of material from the wind by 
the secondary, we assume that, similarly to eq. 5, the mass of 
re-accreted material is proportional to the mass of material 
loss by the white dwarf due to strong wind (M win d)- 



4 RESULTS OF CALCULATIONS 

Since the evolutionary stage of the secondary (hydrogen or 
helium rich) is not observationally determined (more details 
in section 5), we computed four different sequences with hy- 
drogen and helium rich secondaries. Initial parameters of the 
sequences are presented in Table 1. 

Table 1 Initial parameters of computed sequences 



model 


Msg 


M wd 


Pi(RLOF) 


X 


Z 




[M ] 


[M ] 


[d] 






A 


1.4 


0.70 


1.34 


0.7 


0.02 


B 


1.7 


0.85 


1.61 


0.7 


0.02 


C 


1.4 


0.70 


1.26 


0.5 


0.02 


D 


1.7 


0.85 


1.24 


0.5 


0.02 



M cr , 2 = 9.0 x {M wd -0.5), 



(7) 



The sequences C and D present a helium rich SSS chan- 
nel proposed by Hachisu et al. (1999). In Table 2 we present 
results for computed sequences. The orbital parameters are 
given at the moment when orbital period is equal to 1.23 d. 
The effective temperature and luminosity are for subgiant 
star. 

According to our calculations, all models go through 
the short time first nova phase (nl). We have used the grid 
of models calculated by PK95, KP97 and their classifica- 
tion scheme to relate theoretical models to observations of 
nova outburst (for more details see MS98). After that, sys- 
tems enter into the first stable hydrogen burning phase (si), 
which is followed by the wind phase (w). The sequences B, 
C and D have wind phase and second stable hydrogen burn- 
ing phase (s2) . The sequence A does not exhibit wind phase, 
and after stable burning phase this system evolves into sec- 
ond nova phase (n2). After second stable hydrogen burning 
phase sequences B and C evolve through recurrent novae 
phase. In Table 3 we give the duration of each phase. The 
duration of the SSS stage (stable hydrogen burning and wind 
phases) is ranging from several hundred thousand to several 
million years. During the second nova phase, sequence A 
shows behaviour characteristic for slow novae. Sequence D 
avoids second novae phase, and for this system the white 
dwarf mass will exceed Chandrasekhar limit during stable 
hydrogen burning, and a supernova explosion may occur. 

U Sco progenitors can be found with the help of evo- 
lutionary sequences B (hydrogen rich) and C (helium rich). 
After stable hydrogen burning stage both sequences enter 
into recurrent nova phase (Fig.l). In Fig. 2 the evolution of 
the mass accretion rate versus orbital period is shown for 
two sequences B and C. From Fig. 2 we see that both sys- 
tems evolve through the orbital period P or b =1.23 d twice. 



© 0000 RAS, MNRAS 000, 000-000 



4 Ene Ergma et al. 




1.4 1.2 1 

M wd [M ] 

Figure 1. The evolution of the mass transfer rate versus white 
dwarf mass: (a) for hydrogen rich star as donor (sequences A, B 
in Table 1), and (b) for helium rich star as donor (sequences C, 
D in Table 1). The stable hydrogen burning and wind solutions 
are shown. The lines marked 1, 2 and 3 shows lower boundaries 
for recurrent novae, stable hydrogen burning and strong wind 
regions, respectively. The region restricted by lines 1 and 2 shows 
the recurrent nova phase. 



Table 2 Results for computed sequences 



model 


M sg 


M wd 


logT off 


log L/L© 


Msg 




[M ] 


[M ] 


[K] 




[M yr" 1 ] 


A 


0.545 


0.983 


3.669 


0.153 


2.02xl0~ 8 


B 


0.936 


1.304 


3.693 


0.416 


3.58xl0~ 8 


C 


0.603 


1.255 


3.721 


0.388 


4.09xl0~ 8 


D 


1.696 


0.852 


3.957 


1.707 


1.14xl0~ 8 



During the first crossing of P or b = 1.23d line the luminos- 
ity of the secondary is too high (log L/L ~O.9 and 1.35 
for sequences B and C, respectively) and it does not fit the 
observed absolute magnitude My =+3.8 of U Sco. 

For the same two evolutionary sequences, in Fig. 3 we 
present evolution of the orbital period P or b versus mass of 
subgiant M sg (Fig. 3a) and mass of the white dwarf M wd 
(Fig. 3b). From grid of models (PK95 and KP97) we find 
A— 7.6 mag, and P rC c= 23 and 54 yrs for sequences B and 
C, respectively. Therefore, we can conclude that sequence B 
gives the best fit to observing parameters of U Sco. 



S-, 



\ -10 - 



o 
2 



O 




Figure 2. The evolution of the mass transfer rate versus the 
orbital period of the system: (a) for hydrogen rich star as donor 
(sequence B in Table 1,2), and (b) for helium rich star as donor 
(sequence C in Table 1, 2). The arrows show the directions of 
evolution. The vertical (thin dashed) line marks the position of 
the U Sco orbital period. 



Table 3 Time-scales for evolutionary phases 



model 


At 1 

""novae 




^^wind 


At 2 

"Stat 


At 2 

*-* "novae 






[log 


At/yr)] 






A 


4.00 


6.57 






8.17 


B 


5.97 


5.53 


3.70 


6.40 


7.05 


C 


4.76 


5.41 


5.39 


6.17 


7.30 


D 


4.89 


5.02 


5.60 


5.58 





4.1 Chemical composition of the subgiant and 
ejected matter 

Our program is able to follow in detail the evolution of the 
chemical composition of the subgiant and the ejected matter. 
In Table 4 we show isotopic composition in the envelope of 
the subgiant and the ejected matter for two sequences B and 
C, for a moment when the binary system has orbital period 

forf> = l-23 d. 

If we compare the helium rich model C with the hydro- 
gen rich model B we can see that N/C ratio in the envelope 
of the subgiant is higher for helium rich model than hydro- 
gen rich one, but the isotopic ratio 12 C/ 13 C is higher for the 
hydrogen rich model. In the ejected matter both ratios are 
similar. 

Our theoretical calculations show that He/H ratio for 
the matter lost from the system changes during evolution 
from 0.56 to 1.26. For evolutionary sequences B and C and 
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Table 4 Chemical compositions for sequence B and C 



model 


12 C 


13 C 


14 N 


1B N 


16 o 


17 o 


12 C/ 13 C 


N/C 




[xitr 3 ] 


[xlO- 4 ] 


[xlO" 3 ] 


[xlO- 7 ] 


[xlO" 3 ] 


[xlO" 5 ] 






B subgiant 


1.41 


3.05 


3.16 


5.64 


9.77 


0.58 


4.6 


1.8 


B ejecta 


2.81 


11.51 


42.06 


251.60 


1.77 


200.00 


2.4 


10.6 


C subgiant 


0.41 


1.40 


4.53 


2.71 


9.79 


1.17 


2.9 


8.2 


C ejecta 


1.94 


6.97 


37.10 


83.78 


4.94 


68.48 


2.8 


14.0 



V) 

(0 



o 

Oh 



2.5 



.5 



1.5 



.5 



1 1 1 1 1 1 1 1 1 1 1 

L a ) '- 




: b) 

: 1 


i i i i i i i i 


-, , i 





-1 

-2 
-3 



1.5 



M [M D ] 



Figure 3. The evolution of the orbital period of the system as a 
function (a) of the subgiant mass M sg : solid line - sequence C, 
dashed line - sequence B, and (b) of the white dwarf mass M w ^: 
solid line - sequence C, dashed line — sequence B. Horizontal 
thin dashed lines mark orbital period for U Sco, open circle mark 
positions of subgiant and white dwarf masses for sequence B, open 
triangle for sequence C. For upper panel evolution is going from 
left to right, for lower panel from right to left. 
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Figure 4. The evolution of the red giant surface abundances 
of H, He, C, N and O as a function of the subgiant mass M sg : 
upper panel - sequence B, lower panel — sequence C. Vertical thin 
dashed lines mark the position of subgiant mass for sequences B 
(M sg = 0.926 M Q ) and C (M sg = 0.603 M ) respectively (see 
Table 2 for more details). 



for the orbital period of U Sco, this ratio is about 0.7. Our 
He/H determination is well inside observational determina- 
tion which vary between 0.4 and 2 (Anupama & Dewangan 
2000, Williams et al. 1981). 

Figs. 4 and 5 show the evolution of the abundances of 
helium, hydrogen, carbon, nitrogen and oxygen of the sub- 
giant envelope and ejected matter. The vertical thin dashed 
lines show the place where the orbital period is equal to 1.23 
d. In Fig. 5 we also identify the phases of the binary system 
evolution: nl - first short time nova episode, si - first stable 
hydrogen burning phase, w - strong wind phase, s2 - second 
stable hydrogen burning phase and n2 - second nova phase. 

Chemical composition and isotopic analysis may give 
more information about the evolutionary stage of U Sco. 
Possible observational tests are discused in detail by MS98. 



Unfortunately, the subgiant component in U Sco is too faint 
for infrared spectroscopic observations of CO bands in or- 
der to determine the 12 C/ 13 C ratio. However, we think that 
blue domain spectra of U Sco could show some absorption 
structure in the region of 4216A in the CN sequence like 
one observed for DQ Her (Chanan, Nelson & Margon 1978, 
Schneider & Greenstein 1979, Willimas 1983). Analysis of 
this region of the spectra is more complicated because the 
structure of the CH and CN violet system can be affected by 
some absorption features from the disc. However, we suggest 
that since the matter in the accretion disc reflects the chem- 
ical composition of the subgiant star, analysis of the disc 
will give us the information we need if we use observations 
made during the quiescent phase. The chemical analysis of 
the expanding envelope (Anupama & Dewangan 2000) also 
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Figure 5. The evolution of the ejected matter abundances of 
H, He, C, N and O as a function of the subgiant mass M sg : 
upper panel - sequence B, lower panel — sequence C. For more 
explanation see text. 



UV spectrum is predicted to be observable about 10 days 
after the outburst. For hydrogen-rich model (He/H=0.1) 
the supersoft X-ray component is expected to rise till ~ 
50 days after the outburst to a maximum luminosity of ~ 
3xl0 36 ergs _1 (d/kpc) 2 . For helium-rich model (He/H = 2), 
the maximum luminosity is reached about 20 days after the 
optical outburst. Unfortunately, Kahabka et al. (1999) could 
follow the X-ray outburst of U Sco for only 19-20 days af- 
ter outburst. According to Kato (1996), the X-ray luminos- 
ity behaviour depends on the chemical composition. There- 
fore, the evolution of the X-ray luminosity may give the 
important evidence about the chemical composition of the 
accreted matter. 



6 CONCLUSION 

We calculated several evolutionary sequences to reproduce 
orbital and physical parameters of the recurrent novae of 
U Sco-type. We showed that U Sco systems possibly form 
from binaries with a low-mass C-0 white dwarf as accretor. 
Such a system evolves through several observable stages: re- 
current nova, SSS with stable hydrogen burning and SSS 
with strong wind phases. In final phase of evolution a mas- 
sive white dwarf near the Chandrasckar limit is formed. We 
propose that the evolutionary sequence B is able to pro- 
duce binary system with parameters similar to U Sco. Our 
best fitting model has initial parameters: M sg ,i = 1.7 Mq, 
M wdtl = 0.85 M Q and Pi(RLOF) = 1.61d. Based on 
evolutionary model we constracted a detailed theoretical 
model for U Sco. We found that the best fit parameters 



0.94 M Q , M 



are: M a 

M sg = 3.58 x 10~ 8 M Q yr- 1 for P c 



1.31 M , logWLe 
orb = 1.23056 d. 



0.42, 



will give useful information allowing comparison with theo- 
retical models (see Table 4). 



5 DISCUSSION 

Hachisu et al. (1999) proposed a new evolutionary path 
to SNe la, in which the companion star is helium-rich. 
In their model, typical orbital parameters of SNe la pro- 
genitors are: M wd =1.37 M Q , M S9 ~ 1.3 MO, and M sg ~ 
2xl0~ 7 Mq yr _1 . Based on light-curve analysis, Hachisu et 
al. (2000) also constructed a detailed theoretical model for 
U Sco. They found that the best fit parameters are: M ro d~ 
1.37 M Q , and M ag ~ 1.5 M (a range from 0.8 to 2.0 Mq is 
acceptable). 

However, the helium enriched model poses a serious 
problem. Truran et al. (1988) discussed the composition de- 
pendence of thermonuclear runaway models for the recurrent 
novae of U Sco-type. They showed that for M wd = 1.38 M , 
M a g =1.5 x 10~ 8 M yr -1 , L=0.1 L optically bright out- 
bursts are obtained only for matter with He/H<l-2. Above 
results are consistent with our sequences B and C where 
He/H are equal 0.48 and 0.96, respectively. 

There is an alternative explanation of helium enrich- 
ment as observed in U Sco which may also occur due to he- 
lium enriched winds from the white dwarf (Prialnik & Livio 
1995). 

According to Kato (1996), the supersoft component in 
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